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Abstract 

We study self-referential sentences of the type related to the Liar paradox. In particular, we 
consider the problem of assigning consistent fuzzy truth values to collections of self -referential sen- 
tences. We show that the problem can be reduced to the solution of a system of nonlinear equations. 
Furthermore, we prove that, under mild conditions, such a system always has a solution (i.e. a con- 
sistent truth value assignment) and that, for a particular implementation of logical "and", "or" 
and "negation", the "mid-point" solution is always consistent. Next we turn to computational is- 
sues and present several truth-value assignment algorithms; we argue that these algorithms can 
be understood as generalized sequential reasoning. In an Appendix we present a large number of 
examples of self-referential collections (including the Liar and the Strengthened Liar), we formulate 
the corresponding truth value equations and solve them analytically and/ or numerically. 

Keywords. Self-reference, liar paradox, truth, fuzzy logic, nonlinear equations, root finding algo- 
rithms. 

1 Introduction 

Self referential sentences are sentences which talk about themselves. Sentences of this type often 
generate logical paradoxes, the study of which goes back to the ancient Greeks. Many approaches have 
been proposed to neutralize the paradoxes. In this paper we study the problem from the point of fuzzy 
logic. 

The prototypical example of a self-referential sentence is the Liar Sentence: Epimenides the Cretan 
has supposedly uttered the following sentence: 

"All Cretans are liars." 



The Liar sentence generates a paradox. To see this consider whether the sentence is true or false. 
Following a somewhat loose chain of reasoning, let us assume "All Cretans are liars" to mean that 
everything a Cretan says is not true. Since Epimenides is a Cretan, his statement that "All Cretans 
are liars" is not true; if we take "not true" to mean "false", then the opposite of "All Cretans are 
liars" must be true; and if we take this opposite to be "All Cretans are truth-tellers", then what 
Epimenides says must be true, i.e. it is true that "All Cretans are liars" ; but then it must be the case 
that what Epimenides says is false. It seems that we have entered a vicious circle, concluding first that 
Epimenides' statement is false, then that it is true, then again that it is false and so on ad infinitum. 
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The above reasoning is not rigorous, since, for example, "All Cretans are liars" does not necessarily 
mean that all Cretans always utter false statements; similarly the negation of "All Cretans are liars" 
is ^^Some Cretans are truth-tellers" and so on. 

But the paradox also appears under more exact reasoning. Consider the following sentence: 



In this case, if we assume the sentence to be true, then what it says must hold, i.e. the sentence must be 
false. But then its opposite must be true, i.e. it must be true that "This sentence is true" . But then it 
is true that "This sentence is false" and we have again entered an oscillation between two conclusions: 
first that the sentence is false, then that it is true. Notice that we would enter a similar oscillation 
if we started by assuming the sentence to be false; then we would conclude that the sentence is true, 
which would mean that the sentence is false etc. A similar effect can be obtained using two sentences. 
The following pair is the so-called inconsistent dualist: 



If Ai is true, then A2 is also true; but then Ai must be false and so A2 must be false and so on ad 
infinitum. 

Note that self-reference does not necessarily lead to paradox. The following pair is the so-called 
consistent dualist: 



If Ai is true, then A2 is also true, which confirms that Ai is true. In short, accepting that Ai and A2 
are both true is perfectly consistent. However, note that we could equally well assume that Ai is false, 
which would mean that A2 is also false, which confirms that Ai is false. In short, we can also accept 
that Ai and A2 are both false. Is this a problem? We will consider some possible answers in the sequel. 

We have already mentioned that self-referential sentences have been studied extensively, and from 
several different points of view. In this section we will only discuss some work which is directly related 
to the current paper; a more extensive discussion of the literature will be presented in Sectional 

The application of fuzzy logic to the Liar paradox goes back to a paper by Zadeh [SJ|; in summary, 
he resolves the paradox by assigning to the Liar sentence a truth value of 1/2. Following Zadeh's paper, 
several authors have analyzed self-reference using fuzzy logic (for more details see Section ISJ. 

The current paper has been heavily motivated by the work of Grim and his collaborators |in[ I2H 
I22| I41j. Grim considers collections of self-referential sentences and models the fuzzy reasoning process 
as a dynamical system. A method is presented to map each self-referential collection to a dynamical 
system which represents the reasoning process. Each sentence of the self-referential collection has a 
time-evolving fuzzy truth value which corresponds to a state variable of the dynamical system. Grim 
presents several examples of self-referential collections and studies the properties of the corresponding 
dynamical systems. One of the main points of jlUj is that self-referential collections can generate 
oscillating or chaotic dynamical behavior. 

Grim's formulation is an essential starting point for the current paper; but the issues we address 
are rather different. Our main interest is in obtaining consistent truth value assignments, similarly to 
|51j and unlike Grim who concentrates on oscillatory behavior. A more detailed understanding of our 
approach can be obtained by the following outline of the paper. 



This sentence is false. 



Ai 
A2 



'Sentence A2 is true" 
'Sentence Ai is false" . 



Ai 
A2 



'Sentence A2 is true" 
'Sentence Ai is true". 
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In Section [2 we present the logical framework which will be used for the study of self-referential 
sentences. In this we follow very closely Grim's formulation. In Section |31 we reduce the problem of 
consistent truth value assignment to the solution of a system of algebraic^ nonlinear equations. Such a 
system will, in general, possess more than one solution. It is a rather remarkable fact that, for a very 
broad family of self-referential collections, the corresponding equations possess at least one solution, i.e. 
a consistent fuzzy truth value assignment is always possible, under some mild continuity conditions; 
this is the subject of our Proposition |3J In Section H] we turn to the computation of the solutions. 
We consider this as a separate, algorithm- dependent sub-problem. In other words, the same system of 
algebraic equations can be solved by many different algorithms. We consider several such algorithms 
and, returning to Grim's point of view, we study the dynamics of each algorithm. In a certain sense, 
each such algorithm can be understood as a particular reasoning style, and "human-style", sequential 
reasoning of the form "if Ai then A2, if A2 then not-^i etc." is one among many options. In Section 
121 we take a brief look at the literature on self-referential sentences, paradoxes and related topics and 
relate it to our approach. In Section IHl we summarize and discuss our results. Finally, in the Appendix 
we present several specific examples of self-referential collections, formulate the corresponding equations 
and solve them analytically and/or numerically, by using the previously presented algorithms; we also 
compare the behavior of the algorithms. 

2 The Logic Framework 

The general object of our interest is a finite collection of M sentences which talk about each other's 
truth value (i.e. a self-referential collection); in particular we explore the extend to which the internal 
structure of such a system determines the truth values of the sentences. This question will be answered 
mainly in Section [3.21 In this section we introduce the logical framework which is necessary to address 
the problem. This framework is quite similar to the one used in jlUj . 

We start with a finite set Vi, the set of Ist-level elementary sentences (also called variables): 

Vi = {^1,^2,..., .4m}. 

From the 1st level variables we recursively build Si, the set of 1st level sentences (also called logical 
formulas): 

If Am G Vi then Am G Si 
If Si, ^2 G Si then Bi V ^2, A B2, B[ G Si, (1) 

where V, A,' are the logical operators or, and, negation ^. Note that Vi C Si. 
Next we build V2, the set of 2nd level elementary sentences: 

V2 = {"TV(S) = 6" :B£Si,b£ [0,1]} 

where Tr(...) is shorthand for "The truth value of ..."; in other words, "Tr(i?) = 6" means "The truth 
value of B is 6". Finally, we recursively build S2, the set 2nd level sentences: 

If C G V2 then C e S2 
If Di,D2£ S2 then Di V D2, Di A D2, D[ G S2. (2) 

^I.e. static, not time-evolving. 

^We take the informal point of view that precedence of operators, grouping of terms etc. are well understood from 
the context and do not require special explanation. Similalrly, we treat the use of parentheses in a completely informal 
manner. 
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where V, A,' again stand for or, and, negation ^. Note that V2 C 82- 

We will occasionally use the terms elementary truth value assessments for the elements of V2 
and truth value assessments for the elements of 82- Also, we define a special subset of V2 and the 
corresponding subset of 82 as follows. The elementary Boolean truth value assessments are denoted 
by V2 and defined by 

V2 = {"Tr(S) = 6" -.Be Si, be {0,1}}. 
The Boolean truth value assessments are denoted by 82 and defined recursively as follows 

If C7 e V2 then C7 G 83 
If Di,D2e 82 then Di V D2, Di A D2, D[ e §2. 

Again we have V2 C 82- 

Obviously, we can keep building up the hierarchy of sentences, defining V„ in terms of 8„_i and 
Sn in terms of V„; but going up to V2,82 will be sufficient for the purposes of this paper. Let us 
conclude this section by giving some examples of elements from Vi, 81, V2, 82. 

Vi : Ai,A2,...,Am. 

81 : Ai V^3,^2,(^2 A ^4) ... etc. 

V2 : "TV {Ai) = 1", "TV {A't) = 0", "Tr ((^2 V ^3) A ^1) = 0.3" ... etc. 

82 : ["Tr (a;) = 0" A "Tr ((^1 V A^) A A2) = 0.3"] V "Tr (^3) = 0.8" ... etc. 

3 Truth Value Assignment 

Our main goal is to assign truth values to sentences of the 2nd level (i.e. to elements of 82). As it turns 
out we will also need, as an intermediate step, to assign truth values to sentences of the 1st level (i.e. 
to elements of 81). In this section we present two ways of achieving this. Our main interest is in the 
second method, presented in Section r3.21 the so-called implicit truth value assignment, which makes use 
of the self-referential nature of a particular collection. However, for purposes of comparison, in Section 
13.11 we will review the "classical" , explicit method for truth value assignment to non-self-referential 
sentences. 

3.1 Explicit Truth Value Assignment 

The following explicit method of assigning truth values to elements of 81 is well known. We start by 
assigning an arbitrary truth value to every element of Vi (1st level variable). This is equivalent to 
selecting a function x : Vi [0, 1], i.e. £ Vi we have Tr{Am) = x {Am)', for the sake of brevity 

we will henceforth use the simpler notation 

yAm G Vi : Tr (Am) = Xm- 

Next, take any Be 81; it is a logical formula with variables Ai, Am ■ If we replace every occurrence 
of Am with Xm then we obtain a numerical formula containing the variables xi, ■.■,xm and the operators 

■^From the mathematical point of view, the symbols V, A,' denote different operators in Q and in in the first case 
they operate on elements of Si while in the second case they operate on elements of S2; but their logical interpretation 
is the same in both cases. 
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V, a/. To obtain the truth value of B we choose a particular numerical implementation of V, A,' and 
perform the numerical calculations . I.e. 

B = FB{Ai,...,Am) 

Tt{B) = fB{Tr{Ai),...,Tt{AM)) 

or, more concisely, 

Tr (B) = fs {xi, ...,xm) ■ 

Here Fb is a logical formula {Fb ■ Vi — > Si) and /b is the corresponding numerical formula {/b '■ 
[0,1]^ — > [0,1]) obtained by replacing Am with Xm and now understanding the symbols V,A,' as 
numerical operators (specific examples of the procedure appear in the Appendix). The implementation 
of the logical operators V, A,' by numerical operators, namely t-conorms, t-norms and negations, has 
been studied extensively by fuzzy logicians jl7l I28j . Several typical implementations are presented in 
Table 1. 



Family 


X Ay 


X V y 


x' 


Standard 


min(x, y) 


max(x, y) 


1 


— X 


Algebraic 


xy 


X + y — xy 


1 


— X 


Bounded 


max(0, X + y — 


1) 


min(l, X + y) 


1 


— X 




/ X when y = 




/ X when y = \ 






Drastic 


1 y when x = 


'■) 


1 y when x = 1 


1 


— X 




\ else 




V 1 else / 







Table 1 



In the above manner, the truth value assignment originally defined on Vi (i.e. Tv^Am) = Xm, 
m = 1,2, ...,M) has been extended to Si. Now we can use the truth values of 1st level sentences to 
assign truth values to 2nd level elementary sentences as follows. Given a C £ V2, which has the form 

C = "Tr (B) = 6" 

with i3 G Si, 6 G [0, 1], we define 

Tt{C) = 1-\Tt {B)-b\. (3) 

Tr{B) in (jSJ has already been defined, since i? G Si. Note that, according to ©, the maximum truth 
value of C is 1 and it is achieved when Tr{B) = b; the latter is exactly what C says. More generally, the 
truth value of C is a decreasing function of the absolute difference between Tt(B) and b. This certainly 
appears reasonable^. In this manner we can compute the truth value of every C £ V2. Finally, we can 
extend truth values from V2 to S2 in exactly the same manner as we extended truth values from Vi 
to Si. 

Hence starting with a truth value assignment xi,...,xm on Vi we have obtained a truth value 
assignment for every D £ S2, namely 

Tr {D) = fo (xi, ...,xm) 

where Jd ■ [0, 1]*^ ^ [0, 1]. This explicit assignment of truth values does not involve any self-reference 
or circularity: starting with the initial specification of the truth values of Ai, Am ss xi, ...,xm, the 
truth values of all (1st and 2nd level) sentences are obtained as functions of xi, ...,xm- 

*For further justification of ^ see |1U| . Note however, that a number of other functions could be used; a simple 
example is Tr(C) = 1 - (Tr {B) - bf . 
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3.2 Implicit Truth Value Assignment 

Now suppose that we have a collection of M sentences, which talk about the truth values of each other. 
We will show a method by which the information each sentence conveys about the truth value of the 
other sentences (possibly also about its own truth value) can be used to determine (up to a point) the 
truth values in a consistent manner. 

The previously presented framework gives a formal description of sentences which talk about the 
truth values of other sentences. In particular, members of S2 talk about the truth values of members 
of Si 5 Vi. Our approach can be summarized as follows. 

Suppose we are given M self-referential sentences and consider the set S2 which is generated from 
M elementary sentences. We can pick sentences -Di,..., Dm G S2 which have the same structure as 
the original self-referential sentences (examples of the procedure appear in the Appendix). The only 
difference is that the M self-referential sentences talk about each other, while Di,..., Dm talk about 
some elementary, unspecified sentences Ai, Am- However, since Ai, Am are unspecified, we 
can identify A^ with D^ (for m £ {1, ...,M}). Intuitively, this means that D^ says something about 
the truth values of ^1,..., Am, i-e. about the truth values of Di,..., Dm- This is exactly the situation 
which we were trying to model in the first place. While philosophical objections can be raised about 
this type of self-reference, the situation is quite straightforward from the computational point of view 
(as will be seen in the following) and allows the determination of consistent truth values. 

Let us now give the mathematical details. As mentioned at the end of Section 13.11 the truth 
value of every 2nd level sentence D £ S2 (for fixed M and a specific choice of t-norm, t-conorm and 
negation) is a numerical function fo {xi, .-.,xm), the independent variables xi, ---,xm being the truth 
values of Ai, Am ■ To obtain specific truth values by the procedure of Section it is necessary 
to specify xi,-..,xm- Choose a function $ : {1,2, ...,M} S2- <^(1), $(2), <I> (M) are 2nd level 
sentences^ which can also be denoted as Di,..., Dm- Now (for m G {1, ...,M}) Dm is a logical formula 
Fm {Ai, ...,Am)- In other words, we have (for m = 1, 2, M): 

Di=Fi{Au-,AM) 

D2=F2{Ai,...,Am) 

Dm = Fm{Ax,...,Am)- 
Let us form the system of logical equations 

Ai = Di = Fi{Ai,...,Am) 

A2 = D2 = F2{Ai,...,Am) 

(4) 

Am = Dm = Fm {Ai, .-, Am) - 

The "natural" interpretation of Q) is that A^ says (or is, or means) Dm- @ implies that: 

Tr {Ai) = Tr (L»i) = h (Tr {Ai) , TV (Am)) 
Tr (^2) = Tr {D2) = /2 (Tr (^1) , Tr (Am)) 

(5) 

Tr {Am) = Tr (Dm) = Im (TV (Ai) , Tr {Am)) 

^Obviously, in a particular situation we will choose these sentences which have the same structure as the self-referntial 
sentences in which we are interested - see the Appendix. 
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where fm ■ [0, 1] ^ [0, 1] is the numerical formula obtained (by the procedure of Section 13. 1|) from 
Fm- A simpler way to write (jSJ is 

Xl = fi (xi, ...,xm) 
X2 = h (a^i, ■■■■.xm) 

(6) 

XM = Jm {Xl, ■■■,Xm) ■ 

@ is a system of M numerical equations in M unknowns; we will refer to it as the system of truth 
value equations. Note that in general the truth value equations will be nonlinear. 

Depending on the particular <I> used, © may have none, one or more than one solutions in [0, 1]^ . 
Hence, by specifying a particular ^, we obtain a set of possible consistent truth value assignments for 
the 1st level elementary sentences. In other words, every solution of © is a consistent truth value 
assignment. At this point we must consider the possibility that the set of solutions is empty, i.e. that 
there is no consistent truth value assignment. However, as we will soon see, under mild conditions there 
always exists at least one consistent assignment. Assuming that ® has one or more solutions, we can 
choose one of these to assign truth values to the 1st level elementary sentences; next, using exactly the 
same construction as in Section |H3 we can assign truth values to 1st level sentences, then to 2nd level 
elementary sentences and finally to 2nd level sentences. In particular, it is easy to check that at the 
end of the procedure the 2nd level sentences Di, ...,Dm will receive the truth values originally specified 
by the solution of Q - hence the truth value assignment is, indeed, consistent. 

Let us now show that, under mild conditions, every $ function specifies at least one consistent 
truth value assignment. This is the subject of Proposition |21 However, we first need two auxiliary 
propositions. 

Proposition 1 If the implementations o/V,A,' are, respectively, a continuous t-conorm, a continuous 
t-norm and a continuous negation, then /i,/2, ■■■,fM in ® cli^g continuous functions of {xi,X2, ...,xm)- 

Proof. We give a sketch of the proof (we omit the complete proof for the sake of brevity; the basic 
idea is clear enough.). Suppose that V, A,' are a continuous t-conorm, t-norm and negation. Take some 
m G {1,2, ...,M}. Recall that fm{xi,X2, ■■■,xm) =Tr(Z)m)) where Dm S 82- Now, take any B G Si; 
then Tr(i?) = y {xi,X2, ■■■,xm) and y is a finite combination of V,A,' and xi,...,xm, which is clearly a 
continuous function of the vector^ {xi,X2, ...,xm)'^- Furthermore, let C ="Tr(i?) = 6"; then 

Tr (C) = 1 - |Tr (B) - b\ = I - \y (xi, X2, ...,xm) - b\ 

which is also a continuous function of {xi,X2, ...,xm)'^- Since this is true for every B ^ Si and every 
b G [0, 1], we conclude that Tr(C) is a continuous function of x = (xi,X2, ...,xm)'^ for every C G V2. 
Finally, since Dm G S2, and Ti^Dm) is a finite combination of V, A,' and a finite number of terms 
Tr(C7i), Tr(C2), Tr(C7i) (where Ci,C2, -.Cl G V2) it follows that Tr{D m) is a continuous function 

of {xi,X2, ...,XAf)^. ■ 

Proposition 2 Suppose that X is a nonempty, compact, convex set in . If the function f : X ^ X 
is continuous, then there exists at least one fixed point x G X satisfying 

x = f (x) . 

®We use column vectors; the superscript T indicates the transpose. 
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Proof. This the well-known Brouwer's Fixed Point Theorem. Its proof can be found in a number 
of standard texts, for instance in [31 pp. 323-329]. ■ 

Now we can easily prove the existence of consistent truth value assignments. 

Proposition 3 If the implementations o/V, A/ are, respectively, a continuous t-conorm, a continuous 
t-norm and a continuous negation, then (0) has at least one solution x = {xi,X2, ■■■,xm)'^ £ [0, 1]^^. 

Proof. We define the vector function / (xi,X2, ...,xm) as follows: 

f{xi,X2, ...,Xm) = (/l {XI,X2, ...,XA/) ,/2 {xi,X2, .■.,Xm) , ■■; Im {xi,X2, ...,Xm)Y 

where (for m G {1, 2, M}) (xi, X2, ...,xm) is the function appearing in ©. Since fm {xi^X2, ...,xm) 
computes a truth value, we have fm : [0, 1]*^ ^ [0, 1] and hence / : [0, 1]^ — > [0, 1]^ . Furthermore, 
by Proposition n each fm is a continuous function and so / is also a continuous function. Now we can 
apply Proposition 121 with X = [0, 1]*^. ■ 

When we use Boolean truth value assessments, we can prove an additional result about consistent 
truth value assignments. 

Proposition 4 Suppose that in ^ Di,D2, ...,Dm G S2 and the implementations o/V, A,' are, respec- 
tively, max, min and the standard negation. Then ^ admits the solution (1/2, 1/2, 1/2)^. 

Proof. Take any m S {1, 2, M}; then Fm {Ai,A2, Am) is a combination (through V, A,') of a 
finite number of elements Ci, C2, Cl S V2. Take any Ci (with / E {1, 2, L}); it has the form 

Ci = "Tr (Bi) = bi" (7) 

where Bi G S2 and 6; G {0, 1}. The corresponding numerical term will have the form 

Tr{Ci) = l-\Tr{Bi)-bi\. (8) 

or 

zi = l-\yi- bi\ (9) 

where yi =Ti{Bi) and zi =Tr(C/). Now, yi will be a finite combination of xi, xm through max, min 
and negation operators, hence when xi= X2=---= xm= 1/2 we also get yi = 1/2. Then, for 6/ G {0, 1} 
we also get from @ that zi = 1/2. 

Hence every term appearing in fm (1/2, 1/2, 1/2) (the numerical version of Fm {Ai,A2, Am)) 
will be equal to 1/2. Since these terms will be combined with max, min and negation operators it 
follows that fm (1/2, 1/2, 1/2) = 1/2 and this satisfies the m-th truth value equation: 

Xm = \ = fm{l/2,l/2,...,l/2) . (10) 

Since ^ holds for every m G {1, 2, .., M}, it follows that © admits the solution (1/2, 1/2, 1/2)'^. 
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3.3 Discussion 

We have seen that the assignment of consistent truth values to a collection of self-referential propositions 
can be reduced to solving the system of (algebraic) truth value equations. Furthermore, Proposition |31 
shows that, subject to some mild continuity conditions, the truth value equations admit at least one 
solution. 

Proposition|21indicates that the Liar and related self referential collections cease to be paradoxical in 
the context of fuzzy logic. From the mathematical point of view the situation is rather straightforward. 
A system of truth value equations may have no solution in {0, 1}'^; by expanding to [0, 1]*^ the set in 
which we seek solutions we can guarantee that the system always has at least one solution. Of course, 
this result is achieved at the price of admitting fuzzy solutions, i.e. truth values which fall short of 
certainty. 

Next, let us briefly compare explicit and implicit truth value assignment. Explicit assignment is 
straightforward: starting with the initial (and more or less arbitrary) choice of truth values for the 1st 
level variables, the truth values of 1st and 2nd level sentences are uniquely determined. If initial truth 
values were restricted to be either or 1, then the initial choice would essentially be a choice of axioms, 
i.e. the assertion of certain propositions (or their negations); choosing truth values in [0, 1] can be seen 
as a ^^generalized axiomatization" . Generally, in implicit truth value assignment the truth values of 
1st and 2nd level sentences are not uniquely determined; rather a number of possible consistent truth 
values are obtained. There are of course cases (examples will be presented in the Appendix) where 
there is a single consistent truth value assignment. One is tempted to remark that in such cases the 
collection of self-referential sentences is equipped with an implicit axiomatization. 

Let us now return to the system of truth value equations. We have established that it always has 
a solution; but many additional questions can be asked about it. Let us list a few such questions^. 

1. Under what conditions does © possess a unique solution? 

2. Are some solutions "better" than other? For instance, solutions in the interior of the hypercube 
[0, 1]^''^ are "fuzzier" than solutions on the boundary, which in turn are fuzzier than solutions on 
the vertices. What are existence and uniqueness conditions for vertex or boundary solutions? 

3. What is the structure of the set of solutions? For instance: 

(a) are there conditions under which it is a vector space? 

(b) what is its dimension (vector space dimension, Hausdorff dimension etc.)? 

(c) is it equipped with a "natural" order? is it a lattice? 

4. Assuming that some of the above properties are established, are they invariant under different 
implementations of V, A,' and/ or the function Tr(-)? 

One can also ask computationally oriented questions; the most obvious one is: "how to solve the 
truth value equations?" . There is a large number of root finding algorithms that can be used to this 
end and we will discuss some of them in the next section. However, we believe it is useful to distinguish 
between properties of the truth value equations (such as the properties listed above) and algorithm 
properties. We will further discuss this distinction in Section [4.41 ^. 

^We make no attempt to answer these questions in the current paper; they will be the subject of future research. 

*Note that in our formulation the numbers of unknowns and equations are both equal to M. This is a natural 
consequence of the formulation, since both the unknowns and equations are corresponded to the original M self-referential 
sentences. From the mathematical point of view, since every variable is associated to an equation of the form Xm = 
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4 Computational Issues 

We now discuss a number of algorithms which can be used to solve the truth value equations. 



4.1 Root Finding 

Probably the most popular method for solving systems of nonlinear equations is the iterative Newton- 
Raphson algorithm. A description and analysis can be found in many numerical analysis textbooks 
(e.g. HH]). Using 

let us rewrite the truth value equations © as follows. Define (for m G {1,2, ...,M}): 

(2^) — fm (^1 , ) 

and 

h (x) = {hi (x) , /iM (x))'^ . 
Then the truth value equations © are equivalent to 



h (x) = 0. 

Let us define a time varying vector 

X (t) = {xi (t) , ...,XM (i))^ ; 
then the Newton-Raphson method consists in the iteration (t = 0, 1, ...): 

x{t + l)=x{t)- [G (x (t))]-'^ ■ h (x (t)) 



with 



G(x) 



( 



\ 



dhi 


dhi 


8hi 


dxi 


8x2 




dh2 


8h2 


8h2 


dxi 


8x2 


8xM 


dhM 


8hM 


dhM 


8x1 


8x2 


8xM 



J 



Depending on the choice of t-norms and t-conorms, the partial derivatives of hi 



hM may be undefined 

at some points of [0, 1]*^; however at such points G (x) can be approximated by a differentiable function. 

It is usually stated in textbooks that the Newton-Raphson algorithm will converge to a solution x 
if it starts with x (0) "suficiently close to x". But this statement needs some clarifications. 

fm {xi, ...,xm), the number of variables must equal the number of equations. Of course, the m-th function fm will 
generally not depend on all of xi, xm- 

One could also consider more complicated situations, for instance a collection of M sentences which can be separated 
into two groups of sizes AIi and M2(with Mi + M2 = M) having the following property: the sentences of the first group 
talk only about this group, while the sentences of the second group talk about both groups. This is a self-referential 
collection which contains a closed self-referential sub-collection (namely the first group). The sub-collection corresponds 
to a smaller system of truth value equations; the solutions obtained from this sub-system of Mi equations will in general 
be different from the ones obtained from the original system of M equations. 

Another possibility is to consider a collection of A'^ sentences for which truth values are given and fixed and then a 
second collection of M sentences which talk about themselves and the first collection of sentences. This still results 
in a set of M equations which involve M unknown truth values xi,X2,---,xm and A'^ fixed truth values (call them 
parameters) zi, Z2, zjv. Our Proposition |3|still applies and guarantees the existence of a consistent truth value solution 



(x 



1,X2, 



, Xm 



for every choice of zi,Z2, 



, ZN- 
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A fixed point of (jllj) is a point x which has the fohowing property: if for some t we have x (t) = x, 
then for s = 1, 2, ... we have x (t + s) = x (t). For this to hold it is necessary and sufficient that 

[G{x)]-^ ■ h{x) = 0. (12) 

It is clear from (|12j) that every solution x of © is a fixed point of (|llj) but the converse is not necessarily 
true, i.e. ()11() may have fixed points which are not solutions of ©• Furthermore, a fixed point is called 
stable if 

lim X (t) = X. 

t—fOO 

for every choice of x (0) "sufficiently close" to x; otherwise it is called unstable. The iteration of (|11|) 
can, by definition, converge only to stable fixed points, i.e. some solutions of © cannot be obtained 
because they are unstable fixed points. 

There is an additional possibility: Hll|) may converge to a solution of © which lies outside of 
[0, l]^''^. This problem can be addressed by the following modification of Hll|) : 

1. if, for some t and m, (fTT|) gives Xm (t) > 1, then set Xm (t) = 1; 

2. if, for some t and m, (|TT|) gives x„i (t) < 0, then set x„i (t) = 0. 

In conclusion, the Newton-Raphson algorithm may fail to find a consistent truth value assignment 
for several reasons: 

1. the algorithm may fail to converge; 

2. it may converge to a value which does not solve © or lies outside of [0, 1]^; 

3. it may miss a solution which is an unstable fixed point. 

Furthermore, there is no obvious way in which to obtain all the solutions of © using (|1H) ^. 
However, as a practical matter, the numerical experiments of the Appendix indicate that (|11() 
generally converges to a consistent truth value assignment, i.e. a solution of (jH)) belonging to [0, 1]^'^. 

4.2 Inconsistency Minimization 

Root finding can also be formulated (in a standard manner) as a minimization problem. Let us apply 
this approach to truth value assignment. We write the following inconsistency function 

M 

J (Xi, ...,Xm) = ^ {Xm - fm (xi, ...,Xm))^ • (13) 

m=l 

Note that J {xi, xm) is a reasonable measure of the inaccuracy in satisfying the truth value equations. 
For every choice of (xi, xm)"^, J takes a nonnegative value; every solution of © yields a global 
minimum of J. An intuitive interpretation of (fT3|) goes as follows: {xm — fm {xi, xm))'^ is the 
inconsistency of the m-th sentence and X]^=i i^rn — fm {xi, ...,xm))^ is the total inconsistency of the 
self-referential collection. From Proposition |31 we know that there is at least one consistent truth value 
assignment, i.e. a truth value assignment of zero inconsistency. 

^In fact this remark holds not only for Newton-Raphson, but for all numerical root finding algorithms of which we are 
aware. 
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There is a vast number of minimization algorithms; indeed the "standard" sequential reasoning 
of classical Aristotelian logic can be understood as a simplistic inconsistency minimization algorithm 
which attempts to select the optimizing value of one variable at a time. 

In the experiments reported in the Appendix we will use another simple algorithm, namely steepest 
descent minimization. This consists of the following iteration (for t = 1,2, ...): 



The term ^ in 1)14(1 is the gradient of J with respect to x; at points where J is non-differentiable, 
similarly to the Newton-Raphson case, we can use an appropriate approximation by a differentiable 
function. 

For small enough k every step of ()14|) yields a decrease of J. Hence ()14j) will at least converge 
to the neighborhood of a local minimum of J; however this is not enough to guarantee that @ 
is satisfied. Furthermore, there is no guarantee that x (t) as given by ()14|) will always belong to 
[0, 1]^^. However, note that in place of steepest descent one could use a more sophisticated, constrained 
optimization algorithm to minimize some appropriate function of {xi, xm)'^ under the constraints 
of zero inconsistency and staying in [0, 1]^^. 

Let us also note that the "inconsistency minimization" approach is somewhat related to connection- 
ist cognitive models and could conceivably be used by actual human reasoners. However, we are not 
particularly concerned about the "psychological plausibility" of our approach, neither do we propose 
it as a description of actual human reasoning. 

4.3 A Control-Theoretic Algorithm 

Finally, inspired from control theory, we propose a root-finding algorithm which can be summarized 
by the following equation: 



The motivation for (jl5|) can be explained as follows. Choose any m G {1,2, ...,M}. From (|15|) we see 
that the difference Xm {t + 1) — Xm (t) is equal to —k ■ {xm (t) — fm {xi (t) , ...,xm (t)))- The following 
informal convergence argument holds provided that k is small and, consequently, the change between 
Xm {t) and Xm {t + 1) is also small. In that case we can simply write —k-{xm — fm {xi, ...,xm)) and argue 
as follows. If Xm 

(t) is greater than fm {xi (t) , ...,xm {t)), then Xm {t + 1) becomes smaller than Xm (t) 
and (assuming k to be sufficiently small) the difference Xm — fm {xi, ...,xm) gets closer to zero. Similarly, 
if Xm (t) is smaller than fm {xi (t) , ...jXm (t)), then Xm (t + 1) becomes larger than Xm (t) increases and 
the difference Xm — fm {xi, ...,xm) again gets closer to zero. Finally, if {xm — fm {xi, ...,xm)) equals 
zero, then Xm remains unchanged. Hence, unlike Newton-Raphson and Steepest Descent, the fixed 
points of l|15|) are exactly the solutions of the truth value equations. However, it is not guaranteed that 
every fixed point is stable. 

Eq.p5(l is very similar to schemes used for the control of servomechanisms; standard control the- 
oretic methods can be used to investigate the behavior of x{t), i.e. prove rigorously existence of 
stable fixed points, convergence, boundedness etc. (for example, these issues can be investigated by 
constructing a Lyapunov function of ()15() '). 

Similarly to steepest descent, the control theoretic algorithm has a certain degreee of psychological 
plausibility. Namely, we can imagine a human reasoner who adjusts truth values in small increments, 
according to the discrepancy from the values predicted by the truth value equations. 




(14) 




(15) 
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4.4 Dynamical Systems and Reasoning 

As is the case for every iterative algorithm, the three algorithms presented above can be viewed as 
dynamical systems. Hence our approach is somewhat similar to Grim's, since he also has modelled 
reasoning about self-referential sentences as a dynamical system. There is of course an important 
difference: Grim appears to be more interested in oscillatory and chaotic behavior than in convergence 
to consistent truth values. 

While we find the chaotic behavior of Grim's dynamical systems quite interesting, we believe it is 
a rather secondary aspect of self-referential collections. As we have already mentioned, we believe it 
is a property of a particular algorithm rather than of the self-referential collection itself. Of course, 
the algorithms proposed by Grim are somewhat special, in the sense that they are inspired by human 
reasoning. However, the analogy should not be drawn too far. It is rather unlikely that a human, when 
reasoning about self-referential sentences, simulates chaotic dynamical systems with infinite precision 
arithmetic. We do not bring this issue up to criticize Grim's approach^'^. Our point is that both 
Grim's dynamical systems and the algorithms we have presented in Sections I4.1l - H31 can be seen as 
"generalized reasoning systems". In other words, despite our previous remarks about psychological 
plausibility, the connection of such systems to human reasoning is rather slender; however, they may 
reveal interesting aspects of the mathematical structure of self-referential collections. 

5 Bibliographic Remarks 

The literature on self-referential sentences is extensive and originates from many disciplines and points 
of view. Here we only present a small part of this literature: papers which we have found related to 
our investigation. 

There is a vast philosophical literature on self-referential sentences and the related topic of truth. 
This literature goes back to the ancient Greeks. Some pointers to the early literature can be found in 
the books [21 [221 EH 122]. These books also cover recent developments. Regarding the modern literature, 
we must list Tarski's fundamental paper on the truth predicate [Hj. Two other interesting early papers 
are |33[ I34j . Important developments in the 70's were the concept of truth gaps '48' '49' and Kripke's 
theory of truth ^Hl- Some important papers from the 80's and 90's are (StJ-l., 32 . 43^. 

Preliminary concepts of multi-valued logic j36j have been used in connection to the Liar paradox 
already in the Middle Ages. The notion of truth gaps is also closely related to multi-valued logic, as 
is Skyrms' resolution of the Liar j4Uj . The treatment of the Liar by the methods of fuzzy logic was a 
natural development. We have already mentioned Zadeh's paper [HJ which is, as far as we know, the 
first work dealing with the Liar in the context of fuzzy logic. We have also mentioned that Grim's work 
|ini I2H I22| I41j has been a major motivation for the current paper. In the fuzzy logic literature there 
is a considerable amount of work on self-referential sentences and the concept of truth (for example 
|13[ I14j ) . These approaches fall within the mathematical logic tradition and are quite different from 
ours. But some logicians have also addressed computational issues which are related to our concerns. 
For example the concept of fuzzy satisfiability is addressed in |^ [23 EHI ; on the same topic see jl2j . 

We believe that ideas from the areas of nonmonotonic reasoning and belief revision could yield 
fruitful insights on the Liar and related paradoxes, but we are not aware of any work in this direction. 
However the AGU axioms for belief revision (see for example and the book jSj) appear to us very 
relevant to the study of logical paradoxes. Many of these issues are explored in ^21; further interesting 
recent work includes [Sl 1201 • Also, self-reference is explicitly treated in [201 IHOl • 

We consider our approach to self-referential sentences only marginally related to psychology and 

^^Neither has Grim claimed that his model describes actual human reasoning. 
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cognitive science. However the so-called coherence theory (or theories) of truth, a topic at the interface 
between psychology and philosophy, appears related to self-reference. Probably the best known coher- 
ence theory of truth is the one expounded by Thagard in j45j and further elaborated in j4f)| I47j and 
several other publications. Thagard has also applied his theory to reasoning about a set of sentences, 
each pair of which may cohere or incohere. This theory can also be applied to self-referential sentences 
(though we are not aware of work exploring this connection). Thagard has given a computational 
formulation of his theory as a constraint satisfaction problem 46 ; this formulation is rather similar 
to the inconsistency minimization approach we have presented in Section [4.21 The relationship of co- 
herence to belief networks is discussed |47) . It is interesting to note that Thagard has also presented 
a connectionist formulation of the constraint satisfaction problem, which is related to several connec- 
tionist models of nonmonotonic reasoning [26\ I3H l35]. Finally, from our point of view, a particularly 
interesting paper is jSHl) which proposes a fuzzy measure of coherence. The constraint satisfaction 
approach has also been applied in the context of cognitive dissonance 0135]. It would be interesting 
to apply the theory of cognitive dissonance to reasoning about self-referential propositions, but as far 
as we know this has not been done until now. 

Last but not least, the first chapter of Hofstadter's book ^S] has a large collection of self-referential 
sentences. Some of these are of a quite different flavor from the Liar and are not easily formalized; 
however, we recommend the book to the reader for a more general view on the problem of self-reference. 

6 Conclusion 

In this paper we have presented a fuzzy-logic formulation which can describe a large family of collections 
of self-referential sentences; in particular it can accommodate the Liar, the inconsistent dualist, the 
consistent dualist and the strengthened Liar (see the Appendix, Section IA.7|) . In our formulation, 
subject to some mild continuity conditions, every member of this family admits at least one consistent 
assignment of fuzzy truth values. Hence the main contribution of this paper is to expand Zadeh's 
analysis and to thus show that the Liar and related self referential collections cease to be paradoxical 
in the context of fuzzy logic by the simple expedient of expanding the set of possible solutions from 
{0, 1}^^ to [0, 1]^^. If we accept the price of admitting fuzzy solutions, i.e. truth values which fall short 
of certainty, then we can at one stroke resolve a large number of potential paradoxes. 

Furthermore, we have presented several computational approaches to the problem of finding con- 
sistent truth values. In addition to the rather standard approaches of nonlinear equation solving by 
Newton-Raphson and by minimization we have presented an algorithm inspired by control theory 
which, as demonstrated by a number of numerical experiments, appears to strike a good balance be- 
tween speed of execution, convergence and psychological plausibility (see also the Appendix, Section 

Many issues remain open; in the following paragraphs we list several groups of questions which we 
would like to investigate in the future. 

First, there are mathematical questions regarding the solutions of the truth value equations. Exam- 
ples of such questions include existence and uniqueness of boundary and vertex solutions and invariance 
of the solutions under different implementations of the logical operators and the truth function Tr(-). 
Also, there are questions regarding the algorithmic behavior, for instance the stability of a solution 
under a particular algorithm. The complete resolution of these questions appears to be a difhcult prob- 
lem. Perhaps a fruitful first step will be the study of a restricted family of self-referential collections 
(for example the study of self-referential Boolean truth value assessments). 

A second group of questions concerns optimization issues. There is scope for experimentation with 
alternative implementations of the Tr(-) function, alternative inconsistency functions (for example using 
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absolute values rather than squares) and so on. We have also mentioned that truth value assignment 
could be formulated as constrained optimization, i.e. as the minimization of some fuzziness criterion 
subject to zero inconsistency. 

However, we believe that the most interesting variations of the optimization problem are the ones 
which attempt to address the issue of bounded rationality. Whatever the attractions of infinite valued 
logic may be, we find rather implausible that humans reason with a continuum of truth values. A 
(perhaps small) step to increase the relevance of our approach to actual human reasoning would be 
to minimize inconsistency using a finite (and small) number of truth values. This is an intermdeiate 
step between the classical approach of searching for solutions in {0, 1}^ and the fuzzy logic approach 
of searching in [0, 1]^^. While the formulation of the problem is identical to the one we have presented 
here, a major difference is that now there is no guarantee for the existence of a zero inconsistency 
solution. Hence a number of mathematical questions arise, regarding the properties of the minima of 
the inconsistency function; from the computational point of view, the use of a finite set of truth values 
may require the use of entirely different optimization algorithms. 

Let us also note that there are other possibilities for addressing the issue of bounded rationality, 
for example the use of interval- valued fuzzy sets, the description of truth values in terms of tolerance 
classes etc. These issues require further research. 

Finally, we are very interested in the coherence approach to self-referential sentences. Since self- 
referential sentences make claims about each other's truth values it is rather straightforward to setup a 
network with one node per sentence and connections which are either reinforcing or inhbiting (depending 
on what each sentence says about each other). We have performed preliminary experiments based on a 
Boltzmann machine formulation of the problem; we use a variation of Schoch's fuzzy coherence measure 
pS] . These results will be reported in a future publication. 
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A Appendix: Examples and Experiments 

In this appendix we present several examples of our formulation of self-referential collections. Whenever 
possible, we solve the truth value equations analytically; in addition we apply the three algorithms of 
Section 0] to numerically compute truth values. 

A.l Example 1: The Liar 

We present our first example very briefly since it has already been addressed by Zadeh and several 
other authors. The Liar sentence is 

A = "A is false". 
According to our previous remarks, we define 

C = "Tr [A) =0" 

and then set 

A = C 

which implies 

TV {A) = Tr {C) 
Tr(A) = 1 - |Tr(^) - 0| ^ 
X = 1 - |x - 0| 

X = 1 — X (16) 
where we have set x =Tr(A). (|l(i)) is the truth value equation and it has the unique solution x = 1/2. 

A. 2 Example 2: The Inconsistent Dualist 

Our next example is the inconsistent dualist: 

Ai : "^2 is true" (17) 
A2 : "^1 is false". (18) 

In this case the 1st level elementary sentences are Vi = {^1,^2}. We will use two 2nd level 
elementary sentences 

Ci : "TV {A2) = 1" 
C2 : "TV(^i) = 0". 

Note that Ci and C2 are Boolean elementary truth value assessments, i.e. Ci,C2 G V2 C 82- The 
translation of (|T7|) . (fTH|) consists in mapping Ai to Ci and A2 to C2; then we have the logical equations 

Ai = Ci 

A2 = C2 



from which follows 



TV (^1) = TV (Ci) = 1 - |Tr (^2) - 1| 
TV (^2) = TV (C72) = 1 - |Tr (Ai) - 0| 
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Substituting Tr{Am) with Xm {fn = 1,2) we get 

Xi = 1 — \X2 — 1| 
X2 = 1 — |xi — 0| 

and, since xi,X2 G [0, 1], we finally get the truth value equations 



Xi = X2 (19) 
X2 = l-xi. (20) 



Obviously, eqs. ((T ^ - H2U|) have the unique solution x = (xi, 2:2) = (1/2, 1/2). 
The inconsistency function is 



J (Xl, 3:2) = (Xl - X2)^ + (Xl + X2 - 1)^ 



and has partial derivatives 



|^ = ^((x,-X2)V(x, + X2-l)^)=4x,-2 
dJ d ( , n2 / \2\ 

Xl - X2) + (Xl + X2 — 1) I = 4X2 — 2. 



5X2 5X2 

Setting the partial derivatives equal to zero we obtain 

Qj 

= 4x1 - 2 = (21) 

axi 

QJ 

= 4x2 - 2 = (22) 

0x2 

and, obviously, eqs. (PT|) - I|22|) have the unique solution x = (xi, X2) = (1/2, 1/2) which gives the unique 
truth value assignment of zero inconsistency. 

Furthermore, using the steepest descent algorithm we obtain the the dynamical system 

Xl (t + 1) = Xl {t)-k - (4x1 (t) - 2) (23) 
X2 (t + 1) = X2 (t) - A: • (4x2 it) - 2) . (24) 

This is a linear system which can be written in matrix notation as 

X (t + 1) = ^x (i) + h 

with 

/1-4/c \ L_/^2A: 
^"V \-\k ^~\2k 

A has double eigenvalue: Ai = A2 = 1 — 4A:. We have (for i = 1, 2) 

|Ai| = 1 - 4/c < 1 

for sufficiently small k. Hence (|2;-i |) - (|24j) is stable. The unique fixed point x is obtained by solving 

X = Ax + h. 
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The solution is x = (1/2, 1/2). Because of stability we have for every x (0) that 



lim X (t) 



X. 



This result is verified by numerical simulation. Using k = 0.1 and starting with random initial 
conditions we have performed several simulations of ()23() - ()24() . A typical run is illustrated in Figure 
1. Specifically, Fig.l.a illustrates the evolution of the truth values xi (t) (solid line) and X2 (t) (dash- 
dotted line) as computed by the steepest descent minimization algorithm (P5]l -H24 |l with k = 0.1; 
Fig.l.b illustrates the corresponding evolution of inconsistency. 



Fig. l.a 



Fig. l.b 
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We observe that the algorithm gets to the neighborhood of the optimal truth values x in a few steps 
and thereafter oscillates around x; note also that at equilbrium the algorithm yields an inconsistency 
value larger than zero. This behavior is related to the step size k. Using smaller step size we can 
decrease the amplitude of the oscillation and bring the inconsistency closer to zero, at the expense of 
slower convergence; a typical simulation with k = 0.01 appears in Figure 2. Note that inconsistency is 
now a lot closer to zero. 

Fig. 2.a Fig. 2.b 
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Finally, using the control algorithm, we obtain the dynamical system 



xi {t + l)= xi {t)-k- {xi (t) - X2 (t)) (25) 
X2 it + l)= X2 {t)-k- {X2 {t)-l + XI (t)) . (26) 

This is a linear system which can be written in matrix notation as 

x{t + l) = Ax {t) + b 

with 

-^{':.\'-.): ^-{l) 

A has eigenvalues: Ai = —k + 1 + ik, A2 = —k + 1 — ik. We have (for i = 1,2) 

\Xi\ = \/l - 2/c + 2A;2 < 1 
for sufficiently small k. Hence (|25 ]) -(|26 |) is stable. The unique fixed point x is obtained by solving 

X = Ax + b. 

The solution is x = (1/2, 1/2). Because of stability we have for every x (0) that 

lim X (t) = X. 

This result is verified by numerical simulation. Using k = 0.1 and starting with random initial 
conditions we have performed several simulations of ()25() - (|26() . A typical run is illustrated in Figure 
3. Again, Fig. 3. a illustrates the evolution of the truth values xi {t) (solid line) and X2 (t) (dash-dotted 
line) and Fig.3.b illustrates the corresponding inconsistency (as given by H13|) '). 

Fig. 3.a Fig. 3.b 
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We observe that the control algorithm locates the solution of the truth value equations with higher 
accuracy than the steepest descent method and almost as fast. 
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A. 3 Example 3: The Consistent Dualist 

Our next example is the consistent dualist: 

Ai : "^2 is true" 
A2 : Ml is true" . 

Again the 1st level elementary sentences are Vi = {^1,^2}; the 2nd level elementary sentences are 

Ci : "Tr(^2) = 1" 
C2 : "Tr(^i) = 1" 

which are Boolean elementary truth value assessments, i.e. Ci,C2 G V2 Q S2. Setting 

A2=C2 



we obtain 



TV (Ai) = Tv (Ci) = 1 - |Tr {A2) - 1| 
TV {A2) = Tr (C2) = 1 - |Tr (Ai) - 1| 



and then the truth value equations 



which can be written in simple form as 



Xi = 1 — \X2 — 1| 
X2 = 1 — \xi — 1| 



Xi = X2 
X2 = Xi. 



Any vector of the form x = (/?,/?) (/? G [0,1]) is a solution; i.e. there is an infinite number of 
consistent truth value assignments including complete truth (Tr(^i) =Tr(yl2) = 1) and complete 
falsity {Tr:{Ai) =Tic(A2) = 0); in accordance to Proposition 01 (1/2,1/2) is also a solution. 
The inconsistency function is 



and has partial derivatives 



J = {Xi- X2)^ + (xi - X2y 



_ — = 4xi - 4x2 
0x1 

- — = 4x2 - 4x1. 
ax2 



Setting the partial derivatives equal to zero we obtain 

dJ 



^ 4x1 - 4x2 = (27) 
axi 

dJ 

— = 4x2 - 4x1 = (28) 
ax2 



22 



and it is immediate that H27 |) - H28() has the family of solutions x = {P,(3), each of which gives zero 
inconsistency. Furthermore, using the steepest descent algorithm we obtain the linear dynamical system 



xi (t + 1) = xi {t)-k- (4x1 {t) - 4x2 {t)) (29) 
X2 {t + l) = X2 {t)-k- (4x2 {t) - 4x1 {t)) (30) 

which can be written in matrix notation as 

x{t + l) = Ax {t) 

with 

/ 1 - 4/c Ak \ 
\ 4:k l-Ak ) 

A has eigenvalues: Ai = 1, A2 = 1 — 8A;. We have 

|Ai| = 1 and IA2I Vl - S/c < 1 

hence (|29j) ~ (|3()|> will generally be convergent, but will have oscillatory behavior for certain initial con- 
ditions. The fixed points are obtained by solving 

X = Ax. 



The family of solutions is x = (/?, [3). These results are verified by numerical simulation. Using A; = 0.1 
and starting with random initial conditions we have performed several simulations of (|29j) - (|3()|) . Three 
typical runs are illustrated in Figures 4, 5 and 6; the left panels indicate the evolution of the truth 
values and the right panels illustrate the corresponding inconsistency. We observe that the actual 
equilibrium reached depends on the initial conditions. 

Fig. 4. a Fig. 4.b 
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Fig. 5. a 



Fig. 5.b 
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Fig. 6.a 



Fig. 6.b 
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Finally, using the control algorithm, we obtain the dynamical system 

xi (t + 1) = xi {t)-k- (xi {t) - X2 it)) (31) 
X2 (t + 1) = X2 {t)-k- {X2 (t) - Xi (t)) . (32) 

This is a linear system which can be written in matrix notation as 

x{t + l) = Ax (t) 

with 

/ 1 - /c k 
^"V k 1-k 
A has eigenvalues: Ai = l, A2 = l — 2A;. We have for sufficiently small k 

|Ai| = 1 and IA2I Vl - 2k < 1 

hence H31() -(|32 |) can have oscillatory behavior (for certain initial conditions). Also, the fixed points of 
are obtained by solving 

x = Ax. 
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The family of solutions is, as expected, x = /3) . The results are verified by numerical simulation. 
Using A; = 0.1 and starting with random initial conditions we have performed several simulations of 
(|HT|) - (|H^ . Typical runs are illustrated in Figures 7, 8, 9 (again, the left panels indicate the evolution 
of the truth values and the right panels illustrate the corresponding inconsistency). 
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Fig. 7.b 
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Fig. 9.a 



Fig. 9.b 



A. 4 Example 4 

Our third example is somewhat more comphcated but stih involves only V2 sentences: 



Ai 
^3 



"^2 is true and A3 is false " 
"^1 is true and A3 is false " 
"Ai is false" . 



translates to 



A2 



Di 
D3 



C1AC2 
C3AC2 
C4 



where 





"TV {A2) 


= 1" 


C2 


"TV (A3) 


= 0" 


C3 


"TV (Ai) 


= 1" 


C4 


"TV(Ai) 


= 0" 



We will consider two different implemementations of A. 
A. 4.1 A Implemented by Minimum 

If we implement A by the min t-norm, the truth value equations become 

xi = min [x2, (1 - X3)] 
X2 = min [xi, (1 - X3)] 
X3 = 1 - Xl 

and the inconsistency function is 

J = {xi- min [x2, (1 - ^3)])^ + (x2 - min [xi, (1 - ^3)])^ + (X3 - (1 - xi))^ 
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The truth value equations can be solved analytically. From (|38p we obtain 

Xi = l- X3 

and then - (|57j) become 

xi = min [x2, xi] , X2 = min [xi, xi] 

from which follows that 

Xi = X2, X3 = 1 — Xl. 

In other words, the general solution of - is 

x = {P,P,l-p) 

with P G [0)1]- Note that this includes the extremal solutions (1,1,0) and (0,0,1) as well as the 
mid-point solution (1/2, 1/2, 1/2). Also, using ()38() the inconsistency function becomes 

J = (xi - min [x2, xi])^ + (x2 - min [xi, xi])^ + (x3 - (1 - a:i))^ 

which attains the minimum value of for every x = (/?, /3, 1 — /?) with /3 G [0, 1] . For the steepest 
descent algorithm the dynamical system is 



x{t+l) = x{t) - k 



dJ_ 

dx 



(40) 



We cannot write equations for the gradient explicitly, because the expressions min[x2,l — X3] and 
min [xi, 1 — X3] are not everywhere differentiable. However, using an approximate numerical differen- 
tiation we can simulate (|40|) . In Figures 10 and 11 we present results for two typical simulations with 
k = 0.01. 



Fig. lO.a 



Fig. ll.b 
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Fig. ll.a 



Fig. ll.b 
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Let us also present the results of the control theoretic algorithm; the dynamical system is 

xi {t + l) = xi {t) — k ■ {x2 — min [x2 {t) ,1 — X3 (t)]) 
X2 (t + l) = X2 {t) — k- {x2 — min [xi {t),l — X3 (t)]) 

X3 {t + l)= X3 {t)-k- {X3 - (1 - XI (t))) . 

This system is non-linear. In Figures 12 and 13 we present results for two typical simulations with 
A; = 0.1. 
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Fig.l3.a 



Fig.l3.b 
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A.4.2 A Implemented by Product 

Implementing A by the product t-norm we obtain the truth value equations 

Xl = X2 ■ {1 - xs) 

X2 = Xi ■ {1 - Xs) 
X3 = l- Xl 

We can still use xi = 1 — X3 to simplify the truth value equations to 



from which we obtain 
and finally 

This has the solutions 



Xl = X2- Xl, X2 = Xl, X3 = 1 — Xl 



Xl = xl, X2 = xl, X3 = l- Xl 



xi-[l-xl)=0, X2 = x\, X^ = l-Xl. 



(0,0,1), (1,1,0), (-1,1,0); 

the last solution, however, is inadmissible as a truth value assignment. Hence we see that for the same 
self-referential collection, the product implementation of A yields a subset of the solutions obtained 
through the min implementation. 
The inconsistency function is 

J = (xi - a;2 • (1 - xz) f + {x2 - Xl • (1 - 0:3))^ + (0:3 - (1 - Xi) f 

Because the product operator is everywhere differentiable, we can write explicitly the gradient equa- 
tions. We have 



dJ 2 

= —4X1 + 4X2 + 4xiX3 — 4X2X3 + 2x2Xo 

axi 
dJ 

= 6X1 — 4X2 + 2X3 + 4X2X3 — 4xiX3 -I- 2xiXo — 2 

0x2 
dJ 



3X2, 



2x1 + 2x3 + 4xiX2 - 2x1 + 2x1x2, - 2x\ + 2x\x2 - 2. 
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Using the above expressions and the steepest descent algorithm we obtain results of the type presented 
in Figure 14. Note that in the product implementation the steepest descent algorithm requires a 
considerably larger number of steps to converge than for the min implementation. 



Fig. 14.a Fig. 14.b 
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For the the control theoretic algorithm the dynamical system is 

XI (t + 1) = XI {t)-k- {Xi {t) - X2 (t) • (1 - X3 (t))) 
X2 {t+l)= X2 {t)-k- {X2 (t) - XI (t) • (1 - X3 {t))) 
X3 {t + l)= X3 it)-k-{l- X3 {t)) . 

In Figure 15 we present a typical simulation result. Note the fast convergence to equilibirum. 



Fig. 15.a 
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A. 5 Example 5 

Let us now look at an example similar to the previous one, but involving sentences outside of S2: 

Ai : "The truth value of A2 is 0.90 and the truth value of ^3 is 0.20" (41) 
A2 : "The truth value of Ai is 0.80 and the truth value of ^3 is 0.30" (42) 
^3 : "The truth value of Ai is 0.10". (43) 

(jUJ-ljlSl) translates to 

Ai = Di=CiA C2 
A2 = L>2 = C3 A C4 
Ai = Ds = C5 

where 



Ci 


"Tr(^2) 


= 0.90" 


C2 


"TV (^3) 


= 0.20" 


C3 


"Tr(yli) 


= 0.80" 




"TV (^3) 


= 0.30" 


C5 


"Tr(yli) 


= 0.10" 



We will consider two different implemementations of A. 
A. 5.1 A Implemented by Minimum 

If we implement A with the min operator the truth value equations become 

xi = min [1 - 1x2 - 0.90| , 1 - 1x3 - 0.20|] 
X2 = min [1 - \x2 - 0.80| , 1 - |x3 - 0.30|] 
X3 = 1 - |xi - 0.10| . 

These equations cannot be further reduced and while in principle they can be solved analytically 
by distinguishing cases, this requires an inordinate amount of work. Instead, we will solve them 
numerically, using both Newton-Raphson and the control algorithm. Similarly, we will use steepest 
descent to minimize the inconsistency function 

3 

J (Xl, X2, X3) = ^ (Xi, X2, X3) 
m=l 

where 

Ji (xi, X2, X3) = (xi — min [1 — |x2 — 0.90| , 1 — |x3 — 0.20|])^ 
J2 (xi,X2,X3) = (x2 - min [1 - |x2 - 0.80| , 1 - |x3 - 0.30|])^ 
J3 (xi,X2,X3) = (X3 - 1 - |xi - 0.10|)^ ; 



31 



the dynamical system corresponding to the steepest descent algorithm is 

xi {t + l)= xi (t) 

X2 {t + l)= X2 {t) 
X3 {t + l)= X3 [t] 

For a numerical implementation the partial derivatives can be approximated numerically. Finally, for 
the control theoretic algorithm the dynamical system is 

xi {t + l)= xi {t)-k- {xi (t) - min [1 - \x2 (i) - 0.90| , 1 - |x3 {t) - 0.20|]) 
X2 {t + l)= X2 {t)-k- {x2 (t) - min [1 - \xi (i) - 0.90| , 1 - |x3 (i) - 0.20|]) 
X3 {t+l)= X3 {t)-k- {X3 {t)-l + \X3 {t) - 0.201) . 

In Figure 16 we present some simulation results for the Newton-Raphson algorithm, in Figure 
17 we present results for the steepest descent algorithm and in Figure 18 for the control theoretic 
algorithm. Both Newton-Raphson and the control algorithm discover the same solution, namely x = 
(0.95,0.85,0.15)^. Repeated simulations (not presented here) give always the same solution, so it is 
possible that this is the unique consistent truth value assignment for this problem. The steepest descent 
algorithm gets trapped at a local minimum, namely x = (0.56,0.71,0.59)^, which does not yield zero 
inconsistency. This is not accidental; we have noticed after repeated simulations (not presented here) 
that in general the steepest descent algorithm in this case does not yield zero inconsistency, except 
when started quite close to the true solution. Note the very fast convergence of Newton-Raphson: it 
reaches equilibrium in 7 steps, as compared to the approximately 130 steps required by the control 
algorithm. 



-k 



dJ_ 

dxi 
dJ 

dx2 
dJ 
dx3 



Fig. 16.a Fig. 16.b 
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Fig. 17.a 



Fig. 17.b 
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A. 5. 2 A Implemented by Product 

The situation is similar when we implement A by product.. The truth value equations become 

xi = {l- \X2 - 0.90|) • (1 - \X3 - 0.20|) 
X2 = (1 - \X2 - 0.80|) • (1 - |X3 - 0.30|) 
X3 = l- |xi -0.10| 

Wc omit the details which are very similar to previous cases and present the results of numerical 
simulation. 

In Figure 19 we present some simulation results for the Newton-Raphson algorithm, in Figure 20 we 
present results for the steepest descent algorithm and in Figure 21 for the control theoretic algorithm. 
The Newton-Raphson and steepest descent algorithms reach the solution x = (0.6784, 0.7715, 0.4216)'^, 
while the control algorithm reaches the solution x = (0.0473,0.0872,0.9473)"^. Both of these are 
consistent truth value assignments; other simulations (not presented here) have yielded additional 
solutions. Note that in this case the steepest descent algorithm reaches zero inconsistency. Otherwise, 
the remarks previously made for the min implementation also hold for the product implementation. 
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Fig. 19.a 



Fig. 19.b 
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Fig. 20.a 



Fig. 20.b 
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A. 6 Example 6 

Here is a more complicated example: 

Ai : ("Ai has truth value 0.75" and "^2 has truth value 0.35") or "^4 has truth value 1.00' 
A2 : ("^1 or A3 has truth value 1.00") and "A4 has truth value 0.10" 
^3 : "A2 has truth value 0.00" and "A3 has truth value 0.35" 
A4 : "The opposite of Ai has truth value 0.25" 

This translates to 

Ai = (Ci A C2) V C3 
A2 = C4 A C5 

yl3 = Ce A C7 
A4 = Cs 

where 



Ci 


"The 


truth value 


of 


Ai 


is 


0.75" 


C2 


"The 


truth value 


of 


A2 


is 


0.35" 


C3 


"The 


truth value 


of 


Ai 


is 


1.00" 


C4 


"The 


truth value 


of 


Ai 


V A3 is 1.00 


C5 


"The 


truth value 


of 


Ai 


is 


0.10" 




"The 


truth value 


of 


A2 


is 


0.00" 


Cr 


"The 


truth value 


of 


^3 


is 


0.35" 


Cs 


"The 


truth value 


of 


A 


is 


0.25". 



Note that C4 belongs to S2 proper, i.e. it is not an elementary truth value assessment. 

We consider two alternative systems of truth value equations; in the first case we implement V, A by 
the standard implementation; in the second case we implement V, A by the algebraic implementation. 

A. 6.1 A Implemented by Minimum, V Implemented by Maximum 

We present the truth value equations: 

xi = max [min (1 - \xi - 0.75| , 1 - |x2 - 0.35|) , 1 - |x4 - 1.00|] (44) 

X2 = min [1 — |max (xi, X3) — 1.00| , 1 — |x4 — 0.10|] (45) 

X3 = min [1 - |x2 - 0.00| , 1 - |x3 - 0.35|] (46) 

X4 = l-|l-xi-0.25|. (47) 

We omit further details and directly present the results of numerical simulation. 

In Figures 22, 23 we present some simulation results for the Newton-Raphson algorithm, with 
various initial conditions. Note that in both figures the algorithm fails to locate a consistent truth 
value assignment; this is true for most of the simulations we have run, i.e. usually Newton-Raphson 
fails to solve the truth value equations 
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Fig. 22.a 



Fig. 22.b 
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In Figures 24, 25 we present some simulation results for the steepest descent algorithm, with various 
initial conditions. In Figure 24 the algorithm converges to a;=(0. 875, 0.225, 0.675, 0.875)"^ which yields 
zero total inconsistency; this is not the case in Figure 25, where the algorithm converges to a local 
minimum and the truth value equations are not satisfied. This has occurred in several simualtions, i.e. 
the steepest descent algorithm does not reliably solve the truth value equations (|mi - (|47j) . 
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Fig. 24.a 



Fig. 24.b 
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Fig. 25.b 
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In Figures 26, 27 we present some simulation results for the control theoretic algorithm, with various 
initial conditions. The algorithm always converges to x=(0.875, 0.225, 0.675, 0.875)^ which yields zero 
total inconsistency. Indeed, having run a large number of additional simulations we have noticed that 
the control algorithm always converges; furthermore the equilibrium is always the above mentioned x, 
which perhaps indicates that it is the unique solution of the truth value equations (|44j) - (|47|) . 
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Fig. 26.a 



Fig. 26.b 
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A. 6. 2 A Implemented by Product, V Implemented by Extended Sum 

We present the truth value equations: 



Xl = 


(1- 


|xi-0.75|)-(l- 


|X2-0.35|) + (1- 


- |X4- 1.00|) 






(1- 


|xi-0.75|)-(l- 


|X2-0.35|)-(1- 


|X4 - 1.00|) 


(48) 


X2 = 


(1- 


\xi + — X1X3 — 


1.00|) • (1 - |X4 - 


0.101) 


(49) 


3^3 = 


(1- 


|X2-0.00|)-(1- 


|X3- 0.351) 




(50) 


a;4 = 


1 - 


1-xi -0.25| . 






(51) 



We omit further details and directly present the results of numerical simulation. In Figures 28 and 29 
we present some simulation results for the Ncwton-Raphson algorithm, with various initial conditions. 



The algorithm always finds the solution x = (0.9507, 0.2942, 0.5586, 0.7993)-' . In Fi gures 30 and 31 we 

present some simulation results for the steepest descent algorithm, with various initial conditions. Note 
that in Figure 31 the algorithm converges to a local minimum (nonzero) of the total inconsistency. 
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Fig. 28.a 



Fig. 28.b 
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Fig. 31.a 



Fig. 31.b 




Finally, in Figures 32 and 33 we present some simulation results for the control theoretic algorithm, 
with various initial conditions. The algorithm always converges to the following solution of the truth 



value equations (|i8|) -(j5T |l : x 
solution. 



(0.9507, 0.2942, 0.5586, 0.7993) ^ , rendering likely that this is the unique 
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Fig. 33.a 



Fig. 33.b 




A. 7 Example 7: The Strengthened Liar 

This final example does not fall, strictly speaking, within the framework we have presented. However 
it can be treated by a small, straightforward extension of our approach. 

The example is the so-called "Strengthened Liar", a self-referential sentence which has been often 
used as a test of proposed solutions to the Liar paradox [S7|. The Strengthened Liar sentence is 

A = is not true" . (52) 

To treat this and similar sentences in the fuzzy context we must translate it in terms of a membership 
function for the property of being not true. To this end, consider the sentence 

C = "The truth value of A is not a" . (53) 

A possible truth value assignment for (|53jl is 



Tr(C) 



1 when Tr {A) ^ a 
else 



However, this is too strict. Consider the case when a = 1 and Tr(^) = 0.99. Do we really want to 
assign Tr(C) = 1? How about the case Tr(A) = 0.99999? A more reasonable truth value assignment is 

Tr(C7) = |Tr(^) -a| (54) 
which takes the maximum value of 1 when |Tr {A) — a| = 1, i.e. in the cases 

1. Tr(^) = 1 and a = 0; 

2. Tr(yl) = and a = 1. 

Let us accept and set A = C, i.e. 

^="Tr(A)/a". (55) 
(|55|) is more general than (|52|) : to obtain H52() we set a = 1: 

A = "Tr (A) / 1" . 
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Hence, setting x =Tr(^), we must solve the truth value equation 

x = \x — l\ = l — x 

which has the unique solution x = 1/2. 

Obviously this approach can be extended to treat situations which involve statements of the form 

"The truth value of B is not 6" 

for any i? G Si. Hence our framework can be extended defining V2, the set of 2nd level elementary 
sentences to include both sentences of the form 

Tr(B) = 6 

and 

TV {B) ^ h. 

The definition of S2 remains unchanged but, since it depends on the expanded V2, results to an 
expansion of S2 as well. We leave the details to the reader 

A. 8 Discussion of the Algorithms 

Let us close with a short comparison of the three numerical algorithms. We see that the Newton- 
Raphson algorithm converges very quickly but not always to a solution of the truth value equations. 
Regarding steepest descent, its convergence is guaranteed but not always to point of zero inconsistency. 
It appears that the control-theoretic algorithm combines the best properties: practically guaranteed 
convergence to a point of zero inconsistency and in a relatively small number of iterations. Of course, all 
of the above remarks refer to the experiments of the previous sections; further numerical experiments 
(and, perhaps, analytical work) are required to establish their general validity. 
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